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ABSTRACT 


The prediction of torsional buckling of stiffeners is a necessary 






design skill in the planning of any framed structure or vehicle. The 
likelihood of torsional failure is increased when the stiffener is 
Simultaneously loaded from perpendicular directions, a common loading 
condition for the frames of a ship at sea. This thesis examines first 
the case of a single lateral load on I beam and flat bar stiffeners 
and then the instance of combined lateral and axial loadings. An 
energy method is used and formulated are expressions for the critical 


buckling load for a single laterally applied force and a formula 





relating lateral and axial loads and giving an indication of the 


relative weakening of a stiffener caused by simultaneous application 


of both forces. 


Thesis Supervisor: Dr. P.C. Xirouchakis 


Thesis Reader: Dr. J.H. Williams, Jr. 
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NOTATIONS 


Length of stiffener 

Breadth of I beam flange 

Angle of rotation of stiffener in the YZ plane 
Amplitude of cosine function for B 

Height of stiffener shear center above toe 
Warping Constant 

Additional lowering of the load P due to the load 
being applied a distance h/2 above center 

of the web 

Young's Modulus = 2.068 x 10° KN/m 

Shear Modulus = E/2.6 

Height of stiffener 

Moment of inertia of stiffener about web plane 
Local rectangular coordinates at any point x 
on a stiffener 

St. Venant's torsion constant 

Dimensionless parameter 

Moment about the £ axis 

Lateral load 

Critical lateral buckling load 

Axial end load 


Euler buckling load 








Flange thickness 

Web thickness 

Displacements in the x,y, and z directions 
Strain energy 

Work 

Rectangular coordinates originating at 
centroid of stiffener 

Numerical factor 


Membrane strain 


Poisson's Ratıo 
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INTRODUCTION 


| 
As is the case with most large engineering structures, a ship's 


hull and its interior dividing bulkheads and decks are basically 





EEUU ork patterns of stiffeners covered by an outer shell. It is 














principally this lattice of stiffeners which determines the strength 
and structural stability of a vessel. Structural soundness and 
integrity are particularly critical for a ship because this 
engineering structure must operate in what can be the harshest natural 
environment on earth. Very seldom are the frames of a ship's hull 
subject to loading in a single plane; rather, it is most often a 
Be ination of loads caused by the interaction of the various natural 
and man-made forces from different directions which determine the 
total stress on a member. The subject of this thesis is to develop an 
A pression relating the interaction of simultaneous lateral and axial 
Esdosng of a stiffener. 
i 
Two of the more commonly found shipboard stiffeners, the I beam 
and flat bar type, will be examined first under the action of a 
Pentrally applied lateral load and then under the combined actions of 
Ehe central lateral load and a perpindicularly applied end axial load. 
MEE -utrated point loads were chosen because they represent the most 
.imiting and potentially dangerous cases. An energy approach will be 


ised in all instances. First the strain energy of the system will be 


‘ound and then the work done by the loading forces. By conservation 


10 





of energy, the variation of the strain energy will be equated to the 
work done. From the resulting expression will be developed, in the 
case of the single lateral load, a formula for the critical buckling 
load and, in the case of combined loading, a formula relating the 


interaction of the lateral and axial forces. 


It is hoped that study into the problem of stiffener failure 
under more realistic loading conditions such as this mav aid in the 
formulation and upgrading of design codes used in civilian and 
military ship design, and any efforts that could make a ship safer or 
lighter would certainly be worthwhile. To this end it is hoped that 
the ultimate contribution of this thesis 1S a reasonably accurate 


method of predicting, given a value of axial load, how much lateral 


‚load may be sustained by a stiffener without buckling, or vice versa. 





TI 





CHAPTER 1 


DESCRIPTION OF MODELS 


The two types of strengthening members that will be considered 


¡are the I beam type and flat bar type, to be known respectively as 


Model I and Model II. The coordinate system to be employed with both 
| 
models is shown in figure (1-1). 





IN La 


Figure (1-1) Coordinate System 


me origin of the coordinate system will be taken at the center 


of the beam, with u,v, and w being the displacements along the x, y, 


T2 





and z axes. The deflections will be positive in the direction of the 


axes shown in figure (1-1). 


Torsional buckling is that failure mode characterized by a 
deflection and twisting of the stiffener in the YZ plane. This type 


of deformation is shown in figure (1-2). 





Figure (1-2) Characteristics of Torsional Buckling 


1:53 





1.1 Model I 


Model I as stated has been chosen to be a simply supported I 
beam. The case of a beam with symmetric flanges will be considered, 


and a lateral cross section of Model I is shown in figure (1-3). 


A — 


5 


EE AAA LLL v É — 





Figure (1-3) Model I Lateral Cross Section 


The properties of the Model I cross section are listed below: 
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c - Height of shear center above toe (Flange 
thickness t,f considered negligible compared 


to web height h) 


h 
c = 5 
Cw - Warping constant 
Dee) 
te hg 
Cw = 
24 
J - St. Venant's torsion constant 
1 3 3 
= — + 
J 5 (2bt. tw) 
I2 - Moment of inertia about the web plane 


Contribution ot top and bottom 


flanges only) 
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l.2. Model II 


Model II is the flat bar type stiffener as shown in figures (1-4) 


and (1-5). 


\ 


Figure (1-4) Model II - Flat Bar Stiffener 





Figure (1-5) Model II Lateral Cross Section 
| 





The properties of the flat bar cross section are as follows: 


c - Height of shear center above toe 


Cw - Warping Constant 


Cw = 0 
J - St. Venant's torsion constant 
3 
htw 
J = 
3 
fe ioment o: Inertia about Ehe wed plane 
T 
Iz = 
2 
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CHAPTER 2 


Model I - TORSIONAL BUCKLING UNDER POINT LATERAL LOADING 


2.1 Development of the Strain Energy Equation 


nn 


The case of a simply supported I beam of length a subjected to a 


| 


| point lateral load P in the center is illustrated in figure (2-1) 


IN Lx 


Figure (2-1) I Beam Under Point Lateral Load 


The general equation for the strain energy of a beam of length a under 


bending and torsion is given by (refernce l): 


a/2 2 2 > 
D = ee ES T 
Z 2 dx 2 
dx dx 
O 
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Memeidering only that portion of the beam to the right of an arbitrary 
point x, figure (2-2), the forces acting will be a single one P/2 in 
magnitude acting at x = a/2. If a local coordinate system with axes 
jr k, % is assigned to the point x, the force P/2 will produce a 
moment about the X axis equal to: 


P 


M = 5 (a7 2-= > X Y B (2-1) 








P/2 P/2 


Figure (2-2) Model I showing arbitrary point x 


We rOrecesP/2 acting to the right a distance a/2 = x 


Where B again is the angle of twist indicated in figure (1-2). It 
should be noted that B is assumed to be sufficiently sma!l so that 


small angle approximations for sine and cosine are applicable. 





| The expression for 2 = then is (reference l): 
dx 





ILS) 








s . TX a a 
Boi ne Simple choice B = Bo cos 73 will be made. The cosine 
expression is selected because it is assumed that the stiffener will 


be mounted such that it is restrained from twisting at the ends x = 


a/2 and x = -a/2. This is a reasonable assumption in the case of 


stiffeners fixed between decks or frames of a ship. 








2 
i : À d 
Utilyzing the appropriate terms for M, and B, ——— becomes: 
dx 
se E a TX 
= 5 BOO = (2-2) 
dx Z 


Inserting this expression into the formula for strain energy U yields: 


U = a/2 P^Bo* oru TX a x)? it cg (BST) ¢ in? TX 
$ 4EI, a 2 a e a 
O 
+ ECw Bo* ( E E Poe em ax 
a a 
Carrying out the integration it becomes: 
P po? 3, 6 + ne Bor 1,2 a 2 m4 a 
g= ai a + GJ ( ——) " + ECW BO (=) 2 
ET, 1927 d 


2.2 Development of the Virtual Work Equation 


Work is defined as the movement of a force through a distance. 
In the case of Model I then the virtual work done is obtained by 


multiplying the force P times the summation of all the components of 
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motion in the Z direction for elements of the beam between x = 0 and x 
= a/2. It is not necessary to sum vertical components over the entire 
length of the beam, x = -a/2 to x = a/2 , because each half moves the 
same total distance. The expression for the vertical distance moved 
by the centroid of a cross section of the beam at the arbitrary point 


X is (reference 1): 


a/2 





This equation would suffice alone if the force P were applied at the 
centroid of the cross section at X = 0; but since we are considering 
the more realistic case where the load is applied at the outside 
surface of the beam, a distance h/2 above the centroid, the expression 
must be slightly ammended. Due to the twist of the section at X20, 
the load will drop through an additional distance d = h/2(1 - cos ER.) 


as shown in figure (2-3). 
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mate (2-3) Additional lowering of P due to load 


being applied a distance h/2 above centroid 


By the half angle relation, 


T CoS Bio" 2 S rp 


2 (Bro/ 2) 


RR 


p 72 


Therefore d becomes approximately: 


MEA 


=(h/4) BOS SD 


=(h/4) BOS 
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This term is now multiplied by P and added to the work equation to 


give a more accurate expression for the work done. 








a/2 2 
LES pom E O) 
2 2 4 
dx 
O 
a/2 > 
= P Bo cos ( g > E + Ph_ 50° 
2 
ax 4 
O 
| 
| ? ÉS B P a 
| Again we T RI (5— - x) B , and 
| ex Z 


| when applied the work equation becomes: 


2 2 a/2 
| W = M e E — ( > - X ys dx + — aoe 
| 2EI 

Z O 
Integrating this expression gives: 
| 2 2 
W = P Bo a3 ( 6 + > Ph Bo? 
EI, 961 


2.3 Determination of the Critical Buckling Load 


To determine the value of P that will cause failure of the 
detener, the critical buckling load Per, the principle of 


conservation of energy will be employed. This concept states that 
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OE 200!) A 


when a member is gradually loaded, the kinetic energy is zero and the 
work done by external forces is equal to the strain energy (reference 


2). 


To determine ôW and SU calculus of variations with respect to Bo is 
applied to the work and strain energy equations. First to the work 


equation: 








W(Bo + 6Bo) = (Bo + $Bo) * 


23 2 
= W(Bo*) o EE ae 4 E 
En: ee 4 


D 2 
W = 2BoöBo a ( = IR oo 
EI 2 4 
Z 967 


And in lixe manner to the strain energy equation: 








(Da 





2 3 2 
+7 mue EE 
eo) (Bo+sBo) E ¡A un GJ (—) - + ECw(—) - 
EI 2 a 4 a 4 
Z WI ZN 
DA 3 2 2 
SEAT T Tr 
= U (Bo) ^-Bo6Bo Pa Coo e — (GJ + ECw(-) 
EL 2 4a a 
Z 192r 


po 2 2 
i 2 
§U=2Bo6Bo ES T- (GJ * ECw (2) 
ET 2 4a 
Z921 


The two variational expressions may now be equated, and the formula 


for the critical buckling load derived. 
SW = ôU 
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2 3 2 
EXE === ( — + zezh | 
Z 96T 


2 2 2 
Per a 647 T MZ 
2BoSsBo Ex ( > =; at (GJ + ECw (=) 7 





1927 


D. 3 2 2 
Her a 647 Dor T A 10002 _ 
—— = —5! + A E ew GT. | = 0 


This quadratic may be solved for Pcr, and the result expressed in the 








NOE Mm: 
mee = >) -( 2477 E V —2 2 ( 481 E (e E 
TP a GJ 2 Etr? a GJ 
Eo ECw m2 1/2 7 y ET,GJ 
>) MEM 1» = ) ) REC S on (2-4) 
Lm GJ a 


This formula corresponds directly with equation (6-18) of reference 1 


for similar geometry and loading which, in the notation of this 
thesis, 1S: 
y VET GJ 
E R 
2 


a 


Pcr = 


The expression within braces in equation (2-4) equates to the non- 
dimensional factor Y of (2-5). Figure (2-4) shows the close agreement 
between the values of ¥ calculated from equation (2-4) with those 


listed in tabular form in Table 6-5 of reference l. 
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30 


D 
| UNIS ND 
ð from eq (2-4) — —O9- — 
20 a 
a“GJ 
Ew 
| 
10 di) 
| 
\ 
\ 
\ 
9 
a 
ile 
0 10 20 30 40 50 


Ü 


Figure (2-4) Agreement of calculated values of Y from eq (2-4) 


with tabular values listed in reference (1) 
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CHAPTER 3 


MODEL II - TORSIONAL BUCKLING UNDER POINT LATERAL LOADING 


3.1 Development of the Strain Energy Equation 


Model II again was the flat bar type stiffener illustrated in 
figure (1-4). In this chapter the same type of loading situation as 
that of chapter 2, a point load coincident with the Z axis and applied 
at the upper surface of the stiffener, will be examined. A pictorial 
representation of the loading will appear identical to that shown in 


figure (2-1). 


The equation for the strain energy of a flat bar in lateral 
bending and torsion can be obtained in exactly the same manner as that 


for the I beam. The only difference is that the warping rigidity 


term, Cw, vanishes. This is because the equation for Cw of an I beam 
NS. 
2 
teh b 
Cw = 
24 
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and the flat bar stiffener can be likened to an I beam with the widtn 
b of the top and bottom flanges equal to zero. The strain energy 


expression then is: 


a/2 2 2 2 
P Gl aa, dB 
U = \ 2 ( ECT Te ) dx 
dx 
O 








The simple choice for B, Bo cos — , 1S again made and the 
2 
expression for oy is then identical to (2-1): 
dx 
d^v P a TX 
— 2 ull ( 27 x) BO cos E 
dx DT 


The complete strain energy equation is: 


a/2 
U = p cos? zx (2 - x)* + Gs (82% sin” = | dx 


4EI, a 
O 


which when integrated becomes: 


22 2 2 
eee ce s: cB» à 
4EI, 1927 a 4 
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3.2 Development of the Virtual Work Equation 


The work equation is found again by defining the work done as the 
product of the load P times the summation of all vertical components 
of motion over one half of the bar. The expression for the vertical 
distance moved by the centroid at any point along the X axis is 


identical to the I beam case: 





Also, a correction term must once more be added to the work equation 
to account for the load being applied at the upper surface of the 
plate rather than at the centroid. The ultimate result is that the 


work equations for model I and model II are identical: 


2 2 an 
W = SEC CN eos ie (ss = E dx r Ph Bo“ 
DUST a 2 4 
Z 
O 
2 2 2 
EE Oo a> (A e EX po^ 
EI, 96T 4 


e Determination of the Critical Buckling Load 


By setting the variation of the work done equal to the variation 


of the strain energy the value of Pcr is once again determined. Since 
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ee Work expression for the flat plate is identical to that of the I 


beam, so too is the equation for the variation of the work: 


2 3 2 
SW = 2BodéBo Fc. EL + Ph 
El, 96T 4 


While that for ôU is slightly different owing to the abscence of Cw: 


DN 2 
U(Bo + 6Bo) = (Bo + ES E (22) FIGI 2)? :] 
EI, 1921 4 


DN 2 2 
6U = 2BoéBo DON (LE + GJ (=) E 
EI, 1927 


Equating the work and strain energy and solving for Pcr yields: 


SW = SU 
20,3 2 
2BoéBo nouo. E zn 
EI 967 4 
Z 
273 2 2 
T e to eT) ES a 7 
EI 19271 
Z 
De 2 2 
Der ad (© + T) Eun te CI = 0 
EI, 1927 4 4a 


When solved this quadratic gives for Pcr an expression which may De 


written in a form similar to equation (2-2): 
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2 EI u? DIET 
eed 247 h Z » A | 487 h 2 
2 





en) 
6 + ao ao (3-1) 
, a 
For typical stiffeners the quantity > ie several 
2 OTE 
orders of magnitude greater than ( 487 a) (By Eo . Therefore, for 
6 + 12 a GJ 


purposes of a design formula it is sufficiently accurate to ignore 


this latter quantity and rewrite equation (3-1) as: 


P 2 du E ENS 
NE D IPTE oe n | z 





2 6+ 7 6 +7 a GJ a 


(352) 


The formula for Pcr for a flat plate stiffener subjected to a 
point load in the center of the span is given in reference (1) by 


ExUuedon (6-37), which in this thesis' notation is: 








fEI GJ o EI, 
Ber =) 16.94 5 (1 - .87 = ) 
a GJ 


Now to two-decimal accuracy equation (3-2) may be rewritten as: 


es EH. 
Pcr - 17.16 N——E—  (1- .87 


2 
a 





wc 


The formulas are in agreement within 1.33. 





CHAPTER 4 


MODELS I AND II - TORSIONAL BUCKLINS UNDER COMBINED 


LATERAL AND AXIAL POINT LOADING 
Description of Loading Condition 


In chapters 2 and 3 the failure by torsional buckling of a beam 
subject to a lateral load only has been studied; now suppose that an 
additional axial load is added at the ends of the beam as shown in 


figure (4-1). 


A bs 


” 


Figure (4-1) Combined Lateral and Axial Loading of a 


Stiffener 
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| The stiffener is again simply supported and laterally loaded by a 
| 

force P applied at the upper surface in the center of the span, but an 
additional force Q acts at the centroid of each end to compress the 
stiffener inward. The two forces are acting at right angles to each 
bther; P along the Z axis and Q along the X axis. This dual axis type 
idea ing would seem to be much more representative of the kind of 

sing conditions to be encountered by the frames of a ship. For 


"xample, two commonly encountered conditions by a ship at sea in large 


I. are hog and sag, shown in figure (4-2). 


en 


| HOGG ING 





SAGG ING 


i Figure (4-2) Hogging and Sagging Congitions in 


Large Waves 
















Here the ship is supported by the crest of a wave either amidships 
(hogging) or at the ends (sagging). In both cases the vertical 
transverse frames of the ship, a portion of which are shown as dashed 
ines in figure (4-2), will be placed in compression, i.e. a Q force 
applied. At the same time then, the sides of the ship will very 
likely be receiving a lateral force P from other waves in a confused 
sea or in a worst case from collision with another ship, floating 
dbject, or possibly from within the ship from unsecured cargo 
ouffeting about. How do the two forces interact, and what effect does 
the addition of an axial force have upon the capability of the 
stiffener to accept a lateral load? These are the questions this 


thesis will attempt to answer utilizing the strain energy approach 





,» was used in the previous chapters for a lateral load alone. 
| 


a 


| 
ł.2 Development of the Strain Energy Equation 


In solving for the strain energy of a stiffener subjected to both 
Lateral and axial loads the case of model I, an I beam, will be 
considered first; so the energy equation will contain a warping 


megidity term Cw: 
a/2 





2 
pe dB. 2 daB, 2 
) GJ a ) T EC, [——— dx 


) 
dx ase 
O 
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a 
will however, no longer be the same. Taking again 


wain for B, the angle of twist, the choice Bo cos LX is made. The 





expression for 


dx 
an arbitrary point x as shown in figure (2-2), there will still be 


present to the right of x the force P/2 creating a moment about the 


local £ axis equal to the expression in equation (2-1): 


NOSE a a 
Me = > ( 5 XB (2-1) 


There is now though, another force to the right of x adding to this 
moment about? , the axial force O. Looking down upon the stiffener 


from above it is seen that this addition to Mt is equal to the force Q 


times the lever arm -v as shown in figure (4-3): 





Figure(4-3) Additional moment Qv provided by 


axial loading 


ime cotal moment M& then is now: 


P a 
= = Em ES + 
M2 5 (5 x) B Ov 
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ind the expression for 2 E becomes: 


dx 








aoe E M2 
dx“ e 


"CENE UM. a TX. Qv a 
( x ) Bo cos = + EL, (4=1) 





Ve now have a linear second order differential equation which must be 
| 

sOlved for v in order that the strain energy may be found. We are 
Look ing for a particular solution to this differential equation, and 


| 
Enus it will be of the form (reference 3): 


x x MIOS 2 E 4 ` IDS 
cos —— + C., sin + C. x cos -+ C. xXx sin — (4-2) 
a 3 a 4 a 








Differentiating this expression for Vp twice with respect to x yields: 








| 2 
d Vp m lad: TX _ 7 T TX 
5 = (2C, = C,) cos = = (2C4 + = C^) D 
dx 
| 2 2 x 
1 TT AIC: an : m 
- = C4 x cos Eu (3) Cy X sin — (4-3) 


Xf v is a solution to equation (4-1) the following is true: 
p 























2 
d Vp E P c TX Q | 
E ZEI- (5 x) Bo cos = + ET, (C4 eos + C, Su 
TX ) TX 
+ = x COS =~ + C} xsin ) (4-4) 
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By equating coefficients of the trigonometric functions in equations 
(4-3) and (4-4), four equations can be found to be solved for the four 


constants Cy through Cy They are: 


z E er TPBO 
A ad a > DX m 
Z AO EI,T ) 
2 
a PBO p 


O + 2EI T 


Inserting these constants into equations (4-2) and (4-4) gives for E 


ts second derivative: 

















u RES TX m TPBo TX 
Ms = aco — > ) cos Enc ( Z )sin => 
4a Q + 4EI, T (ato x: EI 5^)? 
SS TX 
E a Se ) X cos E (59 
2a Q + 2EI_T 
and 
D 3 
à de "EE 15. 5c ne ERO eos 2S 
2 c ee = D 4a°Q + 4EI 1° 2 
dax Z 
EI TPBO 2 
(Gran + AA ex a 
(aQ+ EL T) 2a Q + 2ET,m 
ES 3 2 E 
Sa O De en 
dx Z 2 4a Q + 4ETT TEV IO) Se 2EI_T 
l nn e 
A B 
3 
er woo e ORS sin T (4-7) 


(a°o + EI, 7)“ 
7 
E 
For simplicity the coefficients in equation (4-7) shall be labeled 


5 





me, and C as indicated, and (4-7) then becomes: 























2 
- D = => (osi (A + Bx) - C sin — 
E 
dx Z 
nd 

acy 2 PBo 2 2 ax 2 202 

— > = ( =—-) °§ cos (MARE DADO IE ) 
EL a 
ax 2 
PBO ; TX 2 2 vc M 
2 EL, (AC + BCx) sin „cos + C sin 5 (4-8) 


Ow with equation (4-8) there is sufficient information to solve the 


train energy equation. 














a/2 2 2 
b.- EI (5-9 0? + c(t « gcw ( 229) ^ | ax 
= dx > ax 
O 
a/2 p^Bo* 2 sex 2 DES 
= —— + x 
\ is COS = (A + 2ABX B x) 
O 
~ 2PBo ( AC + BCX ) sin — Cos 
+ EI c? ne "em DOM sin er 
Z a a 
* ECwBo? ( : )^ cos i = dr (429) 


E~ out the integration of equation (4-9) and reinserting the 
troper quantities for A,B, and C gives as a final expression for 


| 


train energy: 
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eo Oat or Oat E a°o 
E oc Ecco a ox Seen 
Z 1927 DAT 4a" Q + AEI, 7 
6 + x ato 2 
A j 
2 4a Q + 4EI, 7 
p 5 20 ES 
EN M 2 v2 Voce qot 
(4a"Q + AEI, T ) (4a 0 F 4EI T ) 
2022 m 
pP^po^a/ 64EI, T C) T GJ ( Bo 1 
bos 2.4 
(4aQ + AEI, 7 ) 
20 T a 
E ECwBo m (7) (4-10) 


4.3 Development of the Work Equation 


ae work done will again be found from the principle that it is 
force times distance. The work performed by the load P is given by 


the same expresSion used in chapter 2: 





2 
vee \ = Shy (5 ee an - (33) 


d v 
gc 
equation (4-6). Also, there must now also be added the work of the 


Now, however, there will be a new expression for that ob 


Eco at the end x = a/2 as it moves through a distance -u (figure 


4-4). 
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Figure (4-4) Motion of force Q through distance -u 


"The total work done by both loads will be: 
a/2 
E 


EE 








D - x)dx + Ph Bo” + (4-11) 


er E E 4 E 


-All the necessary quantitites in equation (4-11) are known with the 


exception of u . To find this unknown the definition of membrane 


| x=a/2 
Strain (reference 5) will be used: 
du 7 l dw 2 


EIS | 5 (—— ) 


m. za 
E dx 





fe 
DE! 


|Utilyzing inextensibility assumptions this membrane strain is set to 

















Zero: 
u du VR ee) 
EE cu +3 taz? tr a 
Or 
du EN dv 2 1 dw 2 E 
a t5 (ca (4-12) 
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As we are dealing with failure of the stiffener by a torsional 
(twisting) mode augmented by an axial compressive force, it 1s 
expected that deflection in the Z direction will be minor compared 


Therefore, for purposes of a design 





Meee that in the Y direction. 
relationship it will be considered sufficiently accurate to ignore the 


and the relation becomes: 








dw ; 
E term in equation (4-12), 

du Al dv 2 zs 

dx 2? dx ) en 


An expression for u may now be found by solving equation (4-13). Both 


sides are multiplied by dx and then integrated: 





it is necessary to find the 


To find the total work done by O, 
at x=a/2. So 


deflection of the extreme right hand side of the beam, 


the limits of integration must be over the entire length of the beam, 


x=-a/2 to x=a/2; or, twice over half the beam, x=0 to x=a/2: 


a/2 





A 1 dv 2 
pe 725 =? 5 ( dx ) dx (4-14) 


For v equation (4-5) is used, and the first derivative of v with 


respect to x is: 
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av nt ers ESOS DE E F EI EG Noe TX 
dx 4a^Q + 4EI m? 2 eto an 
uc I E 
D E 
a“PBo TX an PBO TX 
+ ( 5 5) cos — 0d 2 wu sco sin = (Ach 

Bang + 2 EI T ASAS MSS 

F G 





Again for simplicity the coefficients in equation (4-15) are labeled 


D,E,F, and G as shown; and equation (4-14) becomes: 




















| a/2 
> NZ TX 2 22 
UE g sin = (D’ - 2DGx + G x ) 
O 
m 2D(F-E) sin — cos + 2G(E-F) x sin cos — 
DNE 2 2 
+ - 
COS ( F 2EF + E) dx (4-16) 
Carrying out the integration of equation (4-16) and replacing 
coefficients gives for the displacement Eaz 
2 
E E 3?p*Bg* ( ON : E ( D 
x=a/2 (4a^Q + AEI 77) i2 4a^Q E AEI v 
8EI ^ > 


R 4a*Q ES 4EI, T 
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All quantities are now known to determine the work done by the 


combined lateral and axial loads. Utilyzing the proper expressions 


























2 
Bor B, ev (eq. 4-6), and u (eq. 4-17), equation (4-11) becomes: 
à 2 x=a/2 
x 
| a/2 
3 TX B a T X 
= P Bo cos z | ZET, ( 5 x ) Bo cos 
O 
3 a^ EI TPBO 
O -a PBo ZO ; TX 
u, |‘, > A Bgm) ae 
i z 4a Q + 4EI T (a Q + EI T ) 
a^PBo TX a Ph 2 
+ ( > 5) SOS S ( 5 sac) ect + Ai Bo 
2a Q + 2EI T 
wo 5 2 ien u^ 
a"P' Bo 18 + 7 ES 28 
(4a^Q + 4EI_ 7°) 12 4a^Q + 4EI «^ 
BEI, 2 
da Q + AEL T 
Integrating equation (4-18) yields for work: 
E e 6 + q? 6 +1 2 6 
EN 1 ec—' DP? 2, 
Z 96 T 24T 4a Q + 4EI, T 
2 5 2 24EI r^ 
P Bo a Q T - 6 2 
BE ——————— ) ( HA pS Mo 
E 4 + | 
mao n AEI T")? a Q AEI T 
2 
SEIT 7 
2 Ph 2 
+ ( ) T 4 BO (4-19) 


4a Q + AEL n“ 


43 








fee Determination of the Critical Buckling Load Equation 


To determine the relationship between the loads P and Q and how 


they will affect buckling of a stiffener, the conservation of energy 


principle will once more be employed. 


First calculus of variations 


with respect to Bo iS applied to equations (4-19) and (4-10) to 


determine as before ¿áW and SU: 











2 2 
SW = 2Bos Bo a — 
Z 961 
p2a°o 
(4a“0 + 4EI T ) 
SEI m 
Z 2 
+ ( q oS 
da Q + 4EI T 
and 
2 3 2 
SU = 2Bo8Bo a [ERUNT - 
Z 1927 
6 ne seg 
E ag E > > 
727 da“Q + 4EI, m 
8a^Q 


- (4 ) | + 


(4a^Q = 4EI 57)" 


id 


2 = 
TT Ir 
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A 5 I ES 5! 
247 4a°Q + 4EI, T 
2 DAET I 
| - 6 Z 





E y 


A) 








12 Neo T 4EI 7f 
Ph 
ER 2 
6 +7 ET 
Ao ©. 
D 4a Q + AEI, T 
| a | SS) 
(4a“0 + 4EI, T ) 
DD 
2 64EI, 7 Q 


nr. o) 
(4a°o + 4EI 7*)^ 


(4-21) 








By conservation of energy 6W = 6U; thus equations (4-20) and (4-21) 
may be equated, and the ultimate relation between P and Q appears as 


the following quadratic with respect to P: 





E un P^a"Q ES. . 
E 1927? (dato + 4EI 07)? 12 
24EI 12 + 8a%0 64(EI 1°)* - 64a“0E1 1° 
om a T r 
da Q + 4EI T (4a“0 + 4EI 7 ) 
Ph a qo 
n a CRE AO (4-22) 


Equation (4-22) presents the relationship between P and Q when they 
are applied to an I beam. If both forces are applied in like manner 


to model II, the flat bar stiffener, the only difference in equation 
T 


4a 


(022) 1s that the final term, ECw( 





3) 1 will be absent since Cw 


meine flat bar is zero. 
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CHAPTER 5 
EVALUATION OF RESULTS OF THE 


BUCKLING LOAD EQUATION 


5.1 Nondimensionalization of the Buckling Load Equation 


It is now our desire to use equation (4-22) to study the 
NEU Go between the lateral and axial loads P and Q and determine 
how changes in the physical dimensions of a stiffener effect its 
ability to withstand against failure. The most interesting comparison 
to be made from equation (4-22) would be a plot of lateral load P 
versus axial load Q, i.e. a means to predict for a given amount of 


lateral load how much simultaneous axial load a stiffener will take 


before buckling and vice versa. 


In order to construct such a plot it will be convenient to first 
uHendimensionalize equation (4-22). An ideal means to 
nondimensionalize P and Q is to express them as a fraction of the 


Euler buckling load (equation 2-5 of reference 1): 


T “1 
Z 


LS oa 
a 
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The Euler critical buckling load Qe is that load applied as in the 


manner shown in figure 5-1 that will just cause in-plane deflection of 


the beam without the twisting effects of torsion. 


ip = 





Figure (5-1) Euler Buckling of a Beam 


Use of the Euler buckling load as a nondimensionalizing parameter is 
particularly well suited to equation (4-22) noticing the freauent 


occurrence in it of the expression 


dao + 4EI_ 1° 


which may be written, if O is taken as positive in tension as: 








-4a*Q + feces = hater ( 1 - ru ) 
The final dimensionless form of (4-22) will be: 
( P y (e E aon (2) 2 
Oe 192 2932 Oe L 
E 37° .Q9 (1)? 2: 5 (9)? (1)? 6 nó Io (1,4 
8 Oe 8 Oe E 4 Qe Ie 


continued next page 
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To three decimal accuracy equation 


mE | En 


Qe 


Q2 


Bos 
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H 


Equation (5-2) 
Borea flat bar 


Germ. 


5.2 Model I - 


the Value 


The resuits of equation (5-2) 


B 


p 2 
d Dm h T J 
(F) + oo a = M ED Es 
) = 0 sl) 
EXC 
E Oe 
(5-1) becomes: 
9. 1,2 1,3 _ 1,4 
z | IE on) 2.467 (7) 
TAS l, 4 P h 
1.234 (7) + 2.467 (p) + 2.467 (57) B 
me ae —— E 
a L, (5-2) 


will be applicable to an I beam, and to make it useful 


stiffener as usual set Cw to zero and delete the final 


Effects of Varving Stiffener Dimensions on 


of the Critical Buckling Load 


for I beams of various sizes are 
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shown graphically in figures (5-2) to (5-5). The points on the graphs 
were derived from computer analysis shown in Appendix A. Along the 
abscissa of each of the graphs is shown the dimensionless ratio of 
axial load Q to the Euler buckling load Qe, while the ordinate gives 
eratio of the critical torsional buckling load Pcr to the Euler 
load. Thus by calculating the Euler buckling load of a stiffener and 
being given a value of axial load Q, we can enter the graphs with a 
value of Q/Qe and obtain a prediction of how much lateral force Pcr, 
applied at the center of the span, that the stiffener will withstand 
before torsional buckling occurs. For example, in figure (5-2) if an 
axial force equal to a value of .4Qe is ea to stiffener 3 with 
h/a ratio equal to .14, then a lateral force of just .16 Qe should 
produce failure of the stiffener by torsional buckling. 

Though no previous design formulas relating simultaneous axial 
and lateral loading of a stiffener to cause torsional buckling were 
found, the results of equation (5-2) do appear reasonable. It is 
possible to verify the predicted results of Pcr/Qe at the intersection 
with the ordinate (Q/Qe=0). For example, again for stiffener 3 of 
figure (5-2) with h/a equal to .14, the complete set of scantlings for 


this I beam are: 


a=3.07m 
n=.43m 
b=.09m 
tf=.02m 


tw=.007m 
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E 


Numbers adjacent to each 
curve correspond to stiffener 
numbers in Appendix A 


Per/Qe 


h/a increasing 





Q/0e 


Figure (5-2) Pcr/Qe vs O/Qe for I beam varying h/a 








h/tw increasing 
Per/Qe BN 


0/0e 


Figure (5-3) Per/0e vs Q/0e for I beam varying h/ tw 


3 








1.0 


Per/Qe 





b/tf increasing 


@7 Oe 


Figure (5-4) Per/Ne vs Q/Qe for I beam varying b/tf 


2 











O 

10 
Eo 

tf/tw increasing 

rer /Qe 12 
EA 
RD. 

0 2 4 6 8 1,0 


0/Qe 


Figure (5-5)  Pcr/Qe vs Q/Oe for I beam varying tf/tw 








Utilyzing these dimensions the warping constant, St. Venant's torsion 


constant, and moment of inertia of the stiffener about the web are: 


2 2 
c —— (oom) oa 107? me 
2 (.09m) (.02m)” + (.43m) (.007m)” SN 
J dd - - 3 - - = 5.292 x 10 m 
2(.09m) 02 EN 
Z r? 


Now the Euler buckling load Qe for this beam calculated from equation 


(2-5) of reference (1) is: 


“EL, 
Ce = EZ 
_ mó(2.068 x 10° KN/m*) (2.430 x 10m”) 
(3.07m) 4 
= 526.2 KN 


To find the critical load for torsional buckling from a lateral load 


alone we use equation (6-18) of reference (1): 


Y V 81,62 
2 
a 


Pcr = 


y is obtained by entering figure (2-4) of this thesis with 
E -7 4 


.à ewm ee 8 KN/m^) (5.292 x 10 'm^) 
ew Fe eig m 


MUN 
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and reading off Y of Sean stnhezseritıcalstrorsional buckling. toad 


er 
per = 13.3 _[(2.068x10 mm^) (2. 430x107 5m) (7.954x10’KN/m2) (5.292x10 7m°) | Y? 
(3.07m) 


peor KN 


A calculated ratio of Pcr/Qe then is: 


Per E 236.1 KN 
Qe 526.2 KN 
= ,449 


The value of Pcr/Qe obtained from equation (5-2) for this same 
stiffener and shown on figure (5-2) for Q/Qe=0 is .455; an agreement 
of 1.3%. Table 5-1 lists the values for Pcr/Qe at 0/0e=0 for all the 
curves of figures (5-2) through (5-5) and shows the good correlation 
between the results of equation (5-2) and the values calculated using 
equations (2-5) and (6-18) of reference (1). All cases are within 2% 
agreement except stiffener 1, figure (5-2), and stiffener 9, figure 
(5-4). In these two cases a problem arises in trying to pick off an 
accurate value of Y from figure (2-4) for use in equation (6-18). We 
are in the vicinity of aa equal to one on figure (2-4). Here the 
W 


values of Y are changing quite rapidly, and a small error in 3 can 


lead to an error of several percentage points in Pcr/Qe. 
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Curve Pcr/Qe at 0/0e=0 Pcr/0e AT 0/0e=0 Percent 


reference (1) Agreement 


equation (5-2) 





Fig 5-2,Beam 1 3836 12.056 6.63 
Beam 2 2596 383 = 

Beam 3 .455 .449 2,3% 

Fig 5-3,Beam 4 .596 2598 . 3% 
Beam 5 2596 2503 25 

Beam 6 S3 259 .5$ 

Fig 5-4,Beam 7 .945 29905 E: 
Beam 8 286 $903 sr 

Beam 9 .486 „sus 5.63 

Fig 5-5,Beam 10 . 748 .734 1.9% 
Beam 11 .586 2593 .5% 

Beam 12 SD -587 2.0% 

Table (5-1) Comparison of values for Pcr/Oe at Q/QezO0 for an I 


beam as calculated bv equation 


equations (2- 


At the right hand side of the curves of figures 


(5-5) the graphs asymptotically approach the x axis, 


Beide = 1. 


(5-2) will have the indeterminate form 1/0. 


5) 


This is because at Q/Qe=1 the expression 1/L in equation 


and 


(6-18) of reference 


(1). 


(5-2) versus values calculated by 


through 


but have no value 


This is to be expected 





o JE —m 


= 


ehouch, because at Q/Oes] the stiffener has failed due to Euler 
buckling, and equation (5-2) can no longer be used to describe the 


relation between Pcr and Q. 


The general trend of all the curves is similar, a rapidly 
decreasing initial slope with the values for Pcr/Qe falling off 
approximateiy 50$ by O/Qe -.37, then tapering off of the slope as the 
curves approach Q/Qe-1. The dimensionless parameters that were varied 
in constructing the graphs were h/a, h/tw, b/tf, and tf/tw. In each 
figure the varied parameter has been doubled twice to construct the 
three curves. It is seen that the ratios which exhibit the greatest 
influence on Pcr/Qe are h/a and b/tf. Varying h/tw has practically no 
effect upon Pcr/Qe, while changing tf/tw at first appears to have a 
fairly significant effect but seems to steady out after tf/tw gets 


above about 2.8. 


5.3 Model II - Effects of Varying Stiffener Dimensions on the 


Value of the Critical Buckling Load 


Figure (5-6) shows the results of equation (5-2) for model II, 
the flat bar stiffener. Once again the values of Pcr/Qe for 0/0e=0 
may be checked against those obtained using the formulas of 
Timoshenko. Consider for example stiffener 13 with h/a=.0523. The 


stiffener will have dimensions and factors as follows: 
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14 


15 


1.4 h/a increasing 


Per/Oe 


Q/Qe 


Figure (5-6) Pcr/Qe vs Q/Qe for Flat Bar varying h/a 





E 





a = 3.059m 
h = .l6m 

tw = .0064m 
Cw = 0 


_ _(.16m) (..0064m) > 8 4 


J 3 = 1.398 x 10 "m 
iom» 0064m)” -9 4 
I = = = 3,495 x 10 m 
zZ 102 


The Euler buckling load from equation (2-5) of reference (1) 13: 


2.068 x 10° KN/m-) (3.495 x 10 m^) 


Qe = > 
(3.059m) 


. 762kKN 


Pcr for the flat bar from equation (6-37) of reference (1) is: 


16.94 VEL GJ h El, 





ga a — (1 - .87 = GF ) 
a 
=1.564 KN 
ker © 1.564 KN E 
Qe u .762 KN = 2.052 


The value of Pcr/Qe obtained from equation (5-2) is 2.079, an 
agreement of 1.3%. The agreement of the values for the other two 


curves are shown in table (5-2). It is seen that there is good 
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agreement in all three cases. 


Curve Pcr/Qe AT O/Qe-0 Pcr/Qe AT Q/Q0e=0 Percent 


equation (5-2) reference (1) Agreement 


FIG 5-6,Beam 13 
FIG 5-6, Beam 14 


FIG 5-6,Beam 15 





TABLE (5-2) Comparison of values of Pcr/Qe at Q/Qes0 for a flat bar 


aS calculated by equation (5-2) versus values calculated by equation 


(6-37) of reference (1). 





As before with the I beam, equation (5-2) is undefined for the 
flat bar at 0/Q0e=1 owing to the quantity 1/L being equal to 1/0; but 
this is acceptable since the bar has undergone Euler buckling at 
Q/Qe=1. The shape and trends of the curves are identical to the I 
beam, and again the curve falls 50$ by about 0/0e=.36. It is 
interesting to note that whereas for the I beam Pcr/Qe is always less 


than one, Pcr/Qe for the flat plate starts out greater than one by a 


60 








factor of approximately two and remains greater than one until 
Q/Oe-.29. This would indicate the much lower relative Qe for the flat 
bar and its much greater susceptibility to failure by Euler buckling 


than the I beam. 


Owing to the geometric simplicity of the flat bar, the only 
dimensionless parameter that may be varied independently is h/a. If 
h/tw is varied while keeping h/a constant, there will be no change in 
equation (5-2). As Cw is zero and the ratio J/Iz for a flat bar is 
always 4, the only term that will alter equation (5-2) is h/a. 


Moreover, as can be seen from figure (5-6), even the changing of h/a 


has very little effect upon Pcr/Qe. 
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CHAPTER 6 


CONCLUSIONS 


In this thesis has been studied the problem of torsional buckling 
of I beam and flat bar type stiffeners under combined lateral and 
axial loads. First examined was the case of lateral loading alone. 

An energy method was employed to determine critical buckling stress 
Per of a stiffener, and equations (2-4) and (3-1) were developed to 
cover lateral loading of I beams and flat bar stiffeners respectively. 
Good agreement was shown between equations (2-4) and (3-1) and the 
published formulas of Timoshenko. One distinct advantage of equation 
(2-4) for an I beam is that it gives an exact formula for the 
coefficient 3 in reference (1) equation (6-18) rather than a tabular 
Iisting of Y for different sizes of stiffeners. As was noted in 
chapter 5, when trying to verify some of the values of Pcr/Qe it 
becomes quite difficult to pick off an accurate value of Y when near 
the bottom of the curve in figure (2-4). In particular, for the two 
cases in Table (5-1) where the percent agreement between Pcr/Qe as 
calculated from this thesis' equations and those of reference (1) were 
6.6% and 5.6%, the differences fall to just .5% and .2% when Y is 


calculated from equation (2-4) rather than read from the graph. 





Next was approached the problem of combined lateral and axial 
loading of a stiffener. Again an energy approach to the problem was 
used, and the ultimate result was equation (4-22) which when 
nondimensionalized using Euler buckling load took the form of equation 
(5-1). Though no other previous formula relating torsional buckling 
under combined loads was found for comparison, the end points of the 
curves of figures (5-2) through (5-6) could be checked and were found 
to be in good agreement with results obtained using other published 
formulas. The trends of all the curves appear reasonable: a constant 
decrease, more rapid at first, of the amount of lateral load a 
stiffener can sustain aS an increasing axial load is applied. As 
another vote of confidence, if O is set to zero in equation (4-22) the 
expression reverts to that of equation (2-4) which was verified 
against equation (6-18) of reference (1). An area for research would 
be to experimentally verify the results of equation (4-22), though it 
would certainly be expensive due to the complexity of a test apparatus 
capable of imposing and measuring Simultaneous lateral and axial 


loads. 


Equation (4-22) is somewhat lengthy for hand calculations, 
Particularly when the constants Cw, H, and Iz must also be found 
before using it. It is not too cumbersome if only one stiffener is 
involved, but if a number of stiffeners are being investigated it does 
lend itself well to a computerized study, especially if an 


investigation is to be done varying stiffener Gimensions. The simple 
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BASIC program with accompanying results shown in Appendix A was 
written and utilyzed to generate the values for constructing figures 


5-2 through 5-6. 
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